Dynamics and rheology of a dilute suspension of vesicles: higher order theory 
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Vesicles under shear flow exhibit various dynamics: tank-treading (tt), tumbling (tb) and 
vacillating-breathing (vb). A consistent higher order theory reveals a direct bifurcation from tt 
to tb if C a = T"7 is small enough (r= vesicle relaxation time towards equilibrium shape, 7=shear 
rate). At larger C a the tb is preceded by the vb mode. For C a > 1 we recover the leading order 
original calculation, where the vb mode coexists with tb. The consistent calculation reveals several 
quantitative discrepancies with recent works, and points to new features. We analyse rheology and 
find that the effective viscosity exhibits a minimum at tt — tb and tt — vb bifurcation points. 
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Vesicles are closed membranes suspended in an aque- 
ous medium. They constitute an interesting starting 
model for the study of dynamics of real cells, such as red 
blood cells. The study of their rheology should capture 
some essential features of blood rheology. 

Under a linear shear flow, a vesicle (where the mem- 
brane is in its fluid state) is known to exhibit a tank- 
treading (tt) motion, while its long axis makes an angle, 
ip < 7r/4, with the flow direction^], In the presence 
of a viscosity contrast A = r)i n /r) (rim and r/ are the in- 
ternal and external viscosities, respectively), ip decreases 
until it vanishes at a critical value of A = A c . For a 
small enough C a = T7 (r is the relaxation time towards 
the equilibrium shape in the absence of an imposed flow, 
7=shear rate) the tt exhibits a saddle-node bifurcation 
towards tumbling (tb) @ . 

Recently, a new type of motion has been predicted 
namely a vacillating breathing (vb) mode: the vesicle's 
long axis undergoes an oscillation (or vacillation) around 
the flow direction, while the shape executes a breath- 
ing motion. Shortly after this theoretical prediction, 
an experimental report on this type of mode has been 
presentedJ5| (trembling denomination was used there) 
and in [6j a qualitatively similar motion called "tran- 
sition motion" in the vicinity of the tt — tb transition 
has been observed. Nevertheless, a detailed experimen- 
tal study of this vb mode would be interesting but has 
not been reported yet. Since then, works providing fur- 
ther understanding^ , or attempting [§, [§] to extend the 
original theory Q to higher order deformation (with the 
aim to account for the experimental observation [5[) have 
been presented. Interesting features have emerg ed@, 
regarding the behavior of the vb mode as a function of 
C a (C a scales out in the leading theory [J]). 

The first aim of this communication is to present the 



result of the consistent theory regarding the higher or- 
der calculation. We find significant differences with a re- 
cent work[j| H[ regarding the form of the evolution equa- 
tion. This implies, in particular, that the location of the 
boundaries separating the various three regimes in pa- 
rameters space is significantly affected. 

A second important report is to investigate how the ef- 
fective viscosity derived recently in[10( is affected by the 
higher order deformation. For a small enough C a the ef- 
fective viscosity of the suspension (as a function of A) still 
exhibits a cusp singularity at the tt — tb bifurcation[Ic| . 
while the cusp becomes a smooth minimum when C a is 
high enough, namely when the tt — vb bifurcation occurs. 

The vesicle suspension is submitted to a linear shear 
flow Vo = (jy, 0,0). The fluid inside and outside the 
vesicle is described by the Stokes equations, with the fol- 
lowing boundary conditions: (i) continuity of the (nor- 
mal and tangential) velocity across the membrane, (ii) 
continuity of normal and tangential stress at the mem- 
brane, (iii) membrane incompressibility. The basic tech- 
nical spirit can be found elsewhere 0, 0, 0, EH Here, 
we merely focus on the results. Lengths are reduced by 
the vesicle radius ro (ro designates the radius of a sphere 
having the same volume). The shape of the vesicle can 
be written in the general case as an infinite series on the 
basis of spherical harmonics y nm 



F nm (t)y nm (8,<p) 
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where e is a small parameter expressing a small deviation 
from a sphere, and 4> are the usual angles in spherical 
coordinates, and F nm (t) is a time dependent amplitude 
(to be determined) of the corresponding spherical har- 
monic. We shall set e = y/A, where A is the membrane 
excess area defined by A = 4ir + A, A being the dimen- 
sionless area of the vesicle. 

Using the expression of spherical harmonics in terms 
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of Cartesian coordinates, r^, we have 
2 

rir fe , (2) 

m=— 2 i^k—x,y,z 

fik are linear combinations of i^m- Since a shear flow 
induces a shape deformation from a sphere which involves 
only second order harmonics (i.e. n = 2), only 3^2m enters 
the calculation^H. 

In recent works[g, Q interesting attempts to deal with 
the higher order deformation (the next step beyond the 
leading order theory [4]) have been made. However, these 
authors take into account only the higher order contribu- 
tion in the Helfrich bending force. We find that there are 
other contributions (especially stemming from the cor- 
responding hydrodynamics response) which are at least 
of the same order. This is the first goal of this commu- 
nication. The next goal is to analyze the far-reaching 
consequences of the present theory. Finally we examine 
some rheological implications. 

The method follows the same strategy as in[4(, and the 
technical details will be reported elsewhere, while we fo- 
cus here on the main results. The extension of the leading 
order calculation jj] to higher order yields, for the ampli- 
tudes of the shape functions fij , the following evolution 
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equation (please compare with Eq.(3) in Ref.[10() 



vt 



20 a 



23A 
4800 



h32 
A- 2 



7 (23A + 32) 
288 (49A + 136) Z 



24 (Z + 6k) 
23A + 32 

— Sd[fipe p j] 



h 



(432A + 1008) k 



7 



(23A + 32) 2 



where e 



, LJ - [diVj + djVil/2 and Sd[b i:ji - 5riJ -r Vji 



3 OijOu\. Zo(t) is the isotropic component of the mem- 
brane tension and is determined by the condition of con- 
stant surface excess area. This constraint provides a 
relation between the f[jS which expresses the fact that 
the shape evolution equations comply with the available 
excess area. Making use of ([3]) fixes then Zq in terms 
of fij (and other parameters, like A, C a ). This rela- 
tion is lengthy and will be listed in an extended paper. 
R = k\/A, where k is the membrane rigidity. Note that 
even though \fK has served formally to make an expan- 
sion of the full equations, it turns out that the amplitudes 
are small enough so that even for A ~ 1 the perturbative 
scheme makes a sense. Our understanding of this fact 
is that A = 1 corresponds to only 8% of relative excess 
area. 

The quantity T>/T>t entering Eq. ^ is the Jaumann 
(or corotational) derivative defined as 
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where M is any second order tensor, D/ Dt is the usual 
material derivative, and w = (Vv — Vv T )/2 is the vortic- 
ity tensor. A fully consistent calculation induces highly 
nonlinear terms. 



Let us first discuss the evolution equation Q and com- 
pare it to recent studies [1, H, Q. For that purpose, it is 
convenient to use another set of variables, namely the ori- 
entation angle tp of the vesicle in the plane of the shear 
and the shape amplitude R by setting F22 = Re~ 21 ^ 0]. 
We may alternatively, instead of R, use the definition 
i?/2v / A = cos(0), as in [8j. We expand the full equa- 
tion in powers of /y and retain terms up to the higher 
(fifth) order in a consistent manner. We then perform a 
straightforward conversion of variables in terms of tp and 
0. We find for and ip the following equations: 



Td t S = -S sin sin 2$ + cos 30 + e A 1 S sin 2$ x 
(cos 40 + cos 20) + e A 2 S sin 2$ cos 20 + . . . . (5) 
S ( cos 2ip 
2 1 cos© 



[l + eA 2 sin 0] - A > + 



where we define, as in 
S 
A 
T 
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and Ai = 



,, V vlWf> ^ = 10/7v/l0A(A- 
2)/(23A + 32). The first term on the right hand side 
of Eq. §5§ corresponds to the leading order theory pre- 
sented in [J], where .... stands for other terms of the 
series. The first term on the right hand side of Eq. 
corresponds to the leading order theory presented in [i[ . 
The first and second term correspond to the situation 
treated in Q where only the higher order contribution 
in the membrane bending force is includedjl, Q. Tak- 
ing the corresponding hydrodynamical response to the 
same order into account (as done here) induce significant 
changes. A new term, for example, is the third term on 
the right hand side of Eq. |(SJ| (proportional to Ai). This 
term is at least of the same order as cos 30. Indeed, the 
term proportional to AiS is of the order of A\C a jt. If 
one has in mind a formal spirit (or a mathematical spirit, 
in that C a is taken of order unity) , then e should be re- 
garded as small. In that case the neglected terms are of 
order 1/e, and are much higher than the retained term 
in the Helfrich energy, namely cos 30 (which is of order 
one). As a natural consequence of this, the so-called sim- 
ilarity equations (put forward in[8J], in that the evolution 
equations contain only 2 independent parameters, S and 
A; while T can be absorbed in a redefinition of time) does 
not hold. Indeed, we have three parameters, which are 
C a , A and A, the excess area (or equivalently e). 

We present now the main outcome following from the 
study of Eq.([3]). We first analyse the tt regime. Figure 
1 presents the orientation angle as a function of A and 
compare the results with previous studies. Instead of a 
square root singularity found for the leading order theory 
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FIG. 1: The angle ipo in the tt regime. Dashed line: the 
leading order theory [j], dotted line: the theory of Ref. §|, the 
full line: the present theory. 
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FIG. 2: The dotted line: theory in Ref.®, the full line: the 
present theory. The same order of discrepancy is found with 
Ref.®. 

(and in the Keller-Skalak regime 13]), the angle crosses 
zero quasi-linearly. A point which is worth of mention 
is that the tt angle becomes negative before the solution 
ceases to exist (signature of the tb regime). Before the 
solution ceases to exist the vb mode takes place, as dis- 
cussed below. 

In [H it was predicted that in the tumbling regime a vb 
mode should take place. This was found to occur as an 
oscillator (like in a conservative system), since the fre- 
quency of oscillation about the fixed point ip = was 
found to be purely imaginary. By including higher order 
terms the frequency acquires a non zero real part [8( , and 
the vb mode becomes a limit cycle (in that all initial con- 
ditions in its domain of existence tend towards a closed 
trajectory in phase space, (ip, a)). As expected from the 
original theory^] the vb mode still occurs in the vicin- 
ity of the tumbling threshold. This happens provided 
that the shape dynamics evolve with time (breathing of 
the shape). C a is a direct measure for the comparison 
between the shape evolution time scale and the shear- 
ing time. The original theoryQ corresponds formally to 
C a = ao, as can be seen from Eq.([5]) and the definition 
of T and S. Including higher order terms leads to the 
appearance of C a in the equation. In what follows, the 
definition of r = Tjajr^/R is adopted. 
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FIG. 3: A snapshot of the vb mode. Solid line: ip = 0. Dotted 
line: ip ~ —0.5. Dashed line ip ~ 0.5 



In Fig. 2 we report on the phase diagram and com- 
pare it to previous theories [H, For small C a we find 
a direct (saddle-node) bifurcation from tt to tb, in agree- 
ment with Ref. [|[ . At C a — > oo we recover the results of 
Ref. [H (in that the vb mode coexists with tb and whether 
one prevails over the other depends on initial conditions; 
this is not shown on the figure 2). At intermediate val- 
ues of C a we find a belt of vb precedingthe tb bifurca- 
tion, in qualitative agreement with [H,I§]. The higher 
order calculation provided here shows significant differ- 
ences with§, as shown Fig. 2. The results presented 
in [1, ^ may be viewed as semi-qualitative given the dis- 
regard of other terms of the same order. Actually, a 
simplistic phenomcnological model captures the main es- 
sential qualitative features of Fig.2[14|. 

Figure 3 shows a snapshot of the vb mode. Note that a 
pure swinging (a terminology usually used for oscillation 
of rigid objects, and adopted in Q) would be impossible 
within the Stokes limit, since this is forbidden by the 
symmetry of the Stokes equation upon time reversal. The 
breathing is a necessary condition for the present mode. 
In the upper half plane (i.e. when ip > 0) the shape 
in the vb regime (dotted line in the Figure) is different 
from the one in the lower plane (full line, rounded shape). 
This asymmetry makes this dynamics possible owing to 
the fact that the two shapes (i.e. for tp > and ip < 0) 
can not be deduced from each other by a simple mirror 
symmetry (albeit this is not easy by necked eyes at first 
look of the figure) with respect to the horizontal axis. 

The basic understanding of the vb mode is as follows. 
First we recall that a shear flow is a sum of a straining 
part along ±7r/4 (which elongates the vesicle for ip > 
and compresses it for ip < 0; see Figure 3) and a rota- 
tional part, tending to make a clockwise tb. Due to the 
membrane fluidity the torque associated with the shear 
is partially transferred to tt of the membrane, so that 
(due to torque balance) the equilibrium angle for tt is 
< ipo < 7t/4. Further, an elongated vesicle tumbles 
more easily than a compressed one[3|. Suppose we are in 
the tt regime (ipo > 0), but in the vicinity of tb, so ipo ~ 
. For small C a the vesicle's response is fast as compared 
to shear, so that its shape is adiabatically slaved to shear 
(a quasi shape-preserving dynamics): a direct bifurcation 
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FIG. 4: The dashed line in the outer graph corresponds to 
the leading order theory where a cusp singularity is observed 
for any C a ■ is the label of the vertical axis is the suspension 
volume fraction. 
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FIG. 5: The first normal stress difference as a function of A 
for two values of C a - 



from tt to tb occurs [3j. When C a ~ 1, the shape does 
not anymore follow adiabatically the shear. When tum- 
bling starts to occur ip becomes slightly negative. There 
the flow compresses the vesicle. Due to this, the applied 
torque is less efficient. The vesicle feels, so to speak, that 
its actual elongation corresponds to the tt regime and not 
to tb. The vesicle returns back to its tt position, where 
ip > 0, and it feels now an elongation (which manifests 
itself on a time scale of the order of 1/7) . Due to elonga- 



tion in this position, tumbling becomes again favorable, 
and the vesicle returns to ip < 0, and so on. We may say 
that the vesicle hesitates or vacillates between tb and tt. 
The compromise is the vb mode. 

Recently a link between the different modes and rhe- 
ology has been presented [l(|- It is thus natural to ask 
how would higher order terms modify the reported pic- 
ture. Eq. ((3|) constitutes a basis for the derivation of the 
constitutive law, as in [l(| • A full study of the derivation 
of the constitutive law will be presented in a future work. 
Here we focus on the effective viscosity as a function of 
A. In the tb and vb regimes we make an average of the 
effective viscosity over a period of oscillation. The results 
are reported on Fig. 4 

We see that at small enough C a , the cusp singularity 
(inset of Fig.4)[T3| persists as in the leading order theory, 
while at larger C a the cusp is smeared out by the fact 
that the transition towards the vb mode does not show a 
singularity as does a saddle-node bifurcation. The same 
reasoning holds for the behavior of N\ (Fig. 5), the normal 
stress difference. For low C a we recover the calculation 
in Ref. 10] (that N\ exhibits a square root singularity, 
while at larger C a the behavior is smoother. 

We have checked that for high enough C a ~ 100 ( 
a quite accessible value in the experiments[6|) the full 
evolution equation produces a coexistence of the vb and tb 
solution, as predicted by the leading order theory (which 
is expected to be recovered formally at high enough C a ). 
We have checked that depending on initial conditions one 
mode prevails over the other (as described inQ). 

Finally, a systematic solution of Eq. ([3]) reveals new in- 
triguing fixed points of the dynamics that do not cor- 
respond to those reported on so far, and they constitute 
presently an important line of inquiry. We hope to report 
along these lines in the future. 

C.V. and G.D. gratefully acknowledges support by the 
European Commission Marie Curie Research Training 
Network MRTN-CT-2004-503661. CM. and T.P. are 
grateful to CNES (Centre National d'Etudes Spatiales) 
and CNRS (ACI "mathematiques de la cellule et du my- 
ocarde") for a financial support. 



[1] M. Kraus, W. Wintz, U. Seifert, and R. Lipowsky, Phys. 

Rev. Lett. 77, 3685 (1996). 
[2] U. Seifert, Eur. Phys. J. B. 8, 405 (1999). 
[3] T. Biben, C. Misbah, Phys. Rev. E 67, 031908 (2003); J. 

Beaucourt, F. Rioual, T. Son, T. Biben and C. Misbah, 

Phys. Rev. E 69, 011906 (2004). 
[4] C. Misbah, Phys. Rev. Lett. 96, 028104 (2006). 
[5] V. Kantsler and V. Steinberg, Phys. Rev. Lett. 96, 

036001 (2006). 

[6] M. Mader, V. Vitkova, M. Abkarian, A. Viallat, T. Pod- 

gorski, Eur. Phys. J. E. 19, 389 (2006). 
[7] P. M. Vlahovska and R. S. Gracia, Phys. Rev. E 75, 

016313 (2007). 



[8] V. V. Lebedev, K. S. T uritsyn, and S. S. Vergeles, 
|arXiv:cond-mat /0702650V l (2007). 

[9] H. Noguchi and G. Gompper, arXiv:cond-mat/0611382 
vl (2006). 

[10] G. Danker and C. Misbah, Phys. Rev. Lett. 98, 088104 
(2007). 

[11] G. Cox, J. Fluid Mech. 37, 601 (1969). 
[12] N.A. Frankel and A. Acrivos J. Fluid Mech. 44 65 (1970). 
[13] S. R. Keller and R. Skallak, J. Fluid Mech. 120, 27 
(1982). 

[14] T. Podgorski (unpublished). 



